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With the full BABAR data sample of 465 x 10 6 BB pairs, we observe the decays B ± -> ipK 1 (l270) ± 
and B ± -> (pK^USO^. We measure the branching fractions (6.1 ± 1.6 ± 1.1) x 10~ 6 and (8.4 ± 
1.8 ± 1.0) x 10~ 6 and the fractions of longitudinal polarization 0.46±£ \l± ° 6 7 and 0.80± ( °f } ± 0.03, 
respectively. We also report on the B — ► (1430)* decay branching fraction of (7.0±1.3±0.9) x 
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10 and several parameters sensitive to CP violation and interference in the above three decays. 
Upper limits are placed on the B ± decay rates to final states with ip and id. (1400) , _K"* (1410) ± , 
K2(1770) ± , or iG(1820) ± . Understanding the observed polarization pattern requires amplitude 
contributions from an uncertain source. 

PACS numbers: 13.25.Hw, 13.88.+e, 11.30.Er 



Measurements of polarization in rare vector-vector B 
meson decay, such as B — > ipK* [2, have revealed 
an unexpectedly large fraction of transverse polariza- 
tion and suggested contributions to the decay amplitude 
which were previously neglected. Decays to other ex- 
cited spin- J kaons Ky can also take place. The differ- 
ential width for a B — > <pKj ' decay has three complex 
amplitudes Aj\, which describe the three helicity states 
A = 0,±1, except when J = 0. The expected hierar- 
chy of the Aj X amplitudes \Aj \ 2 > \Aj+\ 2 > \Aj-\ 2 
is sensitive to the (V — A) structure of the weak inter- 
actions with the left-handed fermion couplings [!, 0, 0] , 
and therefore is sensitive to physics beyond the standard 
model. For example, tensor or scalar interactions would 
violate |Ajo| 2 » |A/+| 2 and the right-handed fermion 
couplings would violate |^4j+| 2 » |^4j-| 2 0- Strong in- 
teraction effects could change these predictions as well, 
but were originally expected to be small 

However, all previous studies have been limited to the 
two-body Kj — > Kn decays, thus considering only the 
spin-parity Kj states with P = (— 1) J . In this paper we 
report the measurement with the three-body final states 
K^p — > Kirn which include P = (— 1) J+1 mesons such as 
K\ and K 2 . We complement these measurements with 
the two-body Kj final states in the £? ± decays and re- 
port polarization in the ^ 1 (1270) ± and tpK^USO)* 1 
final states which have not been seen before. We also 
search for other final states with <p and Kq(1430)^, 
1T 1 (1400) ± , if*(1410) ± , J ftr 2 (1770) ± , or if 2 (1820) ± . 

We use data collected with the BaBar detector 0] at 
the PEP-II e+e" collider. A sample of (465 ± 5) x 10 6 
T(4S) — > BB events was recorded at the the e + e~ center- 
of-mass energy y/s — 10.58 GeV. Momenta of charged 
particles are measured in a tracking system consisting of 
a silicon vertex tracker with five double-sided layers and 
a 40-layer drift chamber, both within the 1.5-T magnetic 
field of a solenoid. Identification of charged particles 
is provided by measurements of the energy loss in the 
tracking devices and by a ring-imaging Chcrenkov detec- 
tor. Photons are detected by a CsI(Tl) electromagnetic 
calorimeter. 

We search for B^ — > ipKj*^ decays using three final 
states of the K { * ]± decay: K ^ , i^Tr , and K^+tt', 
where K s — > 7r + 7r~ and ir° — > 77. We define the two 
helicity angles 0i as the angle between the direction of 
the K or K + meson from K* — > Kir (6>i) or tp — > K + K~ 
(#2) and the direction opposite to the B in the K* or ip 
rest frame. The normal to the three-body decay plane 



for Kj — > i^7T7r is chosen as the analyzer of the Kj 
polarization instead of the direction of K from K*j in the 
two-body decays. We define Hi — cosflj. 

We identify B meson candidates using two kinematic 
variables: m E s = (s/4 - Pb) 1/2 and AE = y/s/2 - E B , 
where (Eb,Pb) is the four-momentum of the B can- 
didate in the e + e~ center-of-mass frame. We require 
m ES > 5.25 GeV and \AE\ < 0.1 GeV (or 0.08 GeV 
for Kj ± — > K ± tt + tt^) GeV. We also require the in- 
variant masses to satisfy 1.1 < vtlk-k < 1.6 GeV, 1.1 < 
tokw < 2.1 GeV, and 0.99 < m K+K - < 1.05 GeV. To 
reject the dominant e + e~ — * light quark-antiquark back- 
ground, we use the angle 9t between the thrust axis of 
the S-candidate decay products and that of the rest of 
the event requiring | cos#t| < 0.8, and a Fisher discrim- 
inant T which combines event-shape parameters Q ■ 

To reduce combinatorial background in the mode 
K*j — > K^ir , we require Hi < 0.6. When more 
than one candidate is reconstructed (7.6% of events with 
Kgir^, 2.9% with K ± tt°, and 14.6% with K^+ir'), we 
select the one whose \ 2 of the charged-track vertex fit 
combined with x 2 of the invariant mass consistency of 
the Kg or n° candidate, is the lowest. We define the b- 
quark flavor sign Q to be opposite to the charge of the 
B meson candidate. 

We use an unbinned extended maximum-likelihood 
fit [ij to extract the event yields rij and the proba- 
bility density function (PDF) parameters, denoted by 
£ and £, to be described below. The index j rep- 
resents the event categories, which include continuum 
background and several B-decay modes. In the S ± — > 
tpKV* 1 — > (K + K~){Kit) topology, the following event 
categories are considered: tpK% (1430) ± , (p(Ktt)q , and 
foiKr)^, where the J p = 0+ (Kir)*^ contribu- 
tion includes both a nonresonant component and the 
K5(1430) ± resonance 0. In the B ± -> <pKj* ]± -> 
(K + K~)(Kmr) topology, we consider t^i^! (1270) ± , 
^(1400)*, ^(1430) ± , ^ J RT*(1410) ± , ^ 2 (1820) ± , 
a nonresonant <pK ± n + TT^ , and /o-HTi(1400) ± contribu- 
tions. In the latter topology, the mode (^^(1770) is 
also considered in place of (,0^2 (1820) ± . In all cases, the 
modes with /o model can account for a possible broad 
non-ip (K + K~) contributions under the tp. 

The extended likelihood is C = exp (— EI A- 

The likelihood Li for candidate i is defined as C{ = 

J2j k n 'j'Pj( Xi '> ^' £)' wnere i s the PDF for variables 
x. = {Hi, H 2 , m Kn ( n ), m K+K -, AE, to es , T, Q}. The 
flavor index k corresponds to the value of Q, that is 
Vj = Vj x SkQ. The C are the polarization parame- 
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ters, only relevant for the signal PDF. The £ parameters 
describe the background or the remaining signal PDFs, 
which are left free to vary in the fit for the combina- 
torial background and are fixed to the values extracted 
from Monte Carlo (MC) simulation [9J and calibration 
B — ► Dtt decays in other cases. 

The signal PDF for a given candidate i is a joint PDF 
for the helicity angles and resonance mass, and the prod- 
uct of the PDFs for each of the remaining variables. 
The helicity part of the signal PDF is the ideal angu- 
lar distribution from Ref. [lCj, multiplied by an empir- 
ical acceptance function Q '(Hi, 'H-i). In the B (pK% 
or ipK2 parameterization, the additional kinematic pa- 
rameters for the decays Kf — > K ± tt + tt~ (such as r\, r%, 
and ro2 in Ref. [I(|) are modeled using the sequential 
two-body decay chains A relativistic spin- J Breit- 
Wigner amplitude parameterization is used for the reso- 



100 



nance masses 



tude p. 



and the J p = 0+ (Ktt)* 



*± 



m 



am- 
plitude is parameterized with the LASS function [8(. The 
nonresonant lpK ± tt + tt~ contribution is modeled through 
sequential if*(892)7r — + Ktttt decay, while the decay 
Kp — ► Ktttt is considered in the systematic uncertainty 
studies. We use a sum of Gaussian functions for the pa- 
rameterization of AE, mES) and T . 

The interference between the J — 2 and (Kn) con- 
tributions is modeled with the term 27Ze(A2oAoo), with 
the three-dimensional angular and mjf, parameteriza- 
tion. We allow an unconstrained flavor-dependent overall 
shift (do + ASq x Q) between the LASS amplitude phase 
and the tensor resonance amplitude phase. The polar- 
ization parameters £ include the fractions of longitudinal 
polarization = |^4jo| 2 /^I^-Ja| 2 in several channels, So, 
and ASo- Similar interference between the ifi(1270) 
and i^i (1400) ± contributions is allowed in the study of 
systematic uncertainties but is not included in the nomi- 
nal fit due to observed dominance of only one mode and 
therefore unconstrained phase of the interference. 

Since the K^IASO)^ meson contributes to all three 
K ^, AT ± 7r°, and K^+tt' final states and (Ktt)*^ 
contributes to two Kir final states in this analysis, we 
consider the total £ as a product of three likelihoods 
constructed for each of the three channels. The corre- 
sponding yields in different channels arc related by the 
relative efficiency. We fit the yields in each charge cat- 
egory k independently and report them in the form of 
the total yield rij — + nj and direct-CP asymmetry 
A C p = (nt -nj)/nj. 

The combinatorial background PDF is the product 
of the PDFs for independent variables and is found to 
describe well both the dominant quark- antiquark back- 
ground and the background from random combinations 
of B tracks. We use polynomials for the PDFs, except 
for mES and T distributions which are parameterized by 
an empirical phase-space function and by Gaussian func- 
tions, respectively. Resonance production occurs in the 
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FIG. 1: Projections onto the variables tties (a), and m K j^ 
(b) for the signal B + — » tp(Kn) and B + — » <p(Kmr) candi- 
dates. Data distributions are shown with a requirement on 
the signal-to-background probability ratio calculated with the 
plotted variable excluded. The solid (dotted) lines show the 
signal-plus-background (combinatorial background) PDF pro- 
jections, while the dashed lines show the full PDF projections 
excluding the signal. 



background and is taken into account in the PDF. 

We observe nonzero B± -> ^X 1 (1270) ± and B ± -> 
( y 9_ft'|(1430) ± yields with significance (excluding sys- 
tematic uncertainties in parentheses) of 5.0(5.3)<7 and 
5.5(6.0)er, respectively. The combined (/3i^i(1270) ± and 
(y5iiri(1400) ± significance is 5.7(6.4)<r. The significance is 
defined as the square root of the change in 2 In £ when 
the yield is constrained to zero in the likelihood C. We 
have tested this procedure with the generated MC sam- 
ples and account for a small observed deviation from the 
one-dimensional \ 2 statistical treatment. 

In Table HI results of the fit are presented, where the 
combined results are obtained from the simultaneous 
fit to the three decay subchannels. In the branching 
fraction calculations we assume K2 — ► i^|(1430)7T and 
B(K*(UW) -> K*tt) = 0.934 ± 0.013 [5]. The signal is 
illustrated in the projection plots in Figs. Q] and [2j where 
in the latter we enhance either the ipKi(1270) ± signal 
(left) or the ( ) o-ft'|(1430) ± signal (right). The nonreso- 
nant K + K~ contribution under the ip is accounted for 
with the B° — > foKi category and its yield 7 ± 16 is con- 
sistent with zero. Similarly, the nonresonant category 
ipKiTTT yield is 148 ± 54 with statistical errors only. 

We vary those parameters in £ not used to model 
combinatorial background within their uncertainties and 
derive the associated systematic errors. Interference 
between the ivi(1270) ± and iYi(1400) ± is one of the 
dominant systematic uncertainties on both yields and is 
modeled with simulated samples. We take the flavor- 
dependent reconstruction efficiency into account in the 
study of asymmetries. The biases from the finite resolu- 
tion of the angle measurement, the dilution due to the 
presence of false combinations, and other imperfections 
in the signal PDF model are estimated with MC simula- 
tion. Additional systematic uncertainty originates from 
possible B background, where we estimate from MC sim- 
ulation that only a few events can fall in the signal region. 

The (/?i4T2(1770) ± yield is not considered in the nom- 
inal fit due to large correlation with (1820) ± . But 
we substitute the ^(1820) resonance for the ^(1770) 
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TABLE I: Results: the reconstruction efficiency e rC co; the total efficiency e, including the daughter branching fractions [jj; 
the number of signal events n s i g ; significance S; fraction of longitudinal polarization fi; the branching fraction B; and the 
flavor asymmetry Acp- The branching fraction B(B ± — > tp^Kir)^ ) refers to the coherent sum \A ICS ± yl n on-rcs| 2 of resonant 
and nonresonant J p = + Kit components || and is quoted for rriKn < 1-6 GeV, while the B{B il — > ip/fo (1430^) is derived 
from it by integrating separately the Breit-Wigner formula of the resonant |^4 rcs | 2 Kn component [j§] without ttiktv restriction. 
When several subchannels contribute, yield and efficiency are quoted for each subchannel. The 90% confidence level upper 
limit on B is quoted with the central values and errors in parentheses. The insert shows two interference parameters So and 
A5q for <pA'2(1430) ± and ^(Kn)^. The tpA'2(1770) ± yield is not considered in the nominal fit and the value indicated with f 
is obtained with these cpA'2(1820) ± yield constrained to zero. The systematic errors are quoted last. 



Mode 


Srcco (%) 


s(%) 




n s i g (events) 


S(a) 




h 


B (10" 6 ) Acp 




c^)A'i(1270) ± 


25.4 ± 1.4 


4.07 ±0 


51 


116 ±26 t\l 


5.0 


n ac+0.12 +0.06 
V-<iV_ 013 _ .07 


6.1 ±1.6 ±1.1 ±0.15 ±0.19 ±0.05 


^(MOO)* 


24.6 ± 1.3 


5.19 ±0 


44 


7 ±39 ±18 


0.2 






<3.2 (0.3 ±1.6 ±0.7) 




<^A'2*(1430) ± 


11.9 ±0.6 


3.34 ±0 
0.64 ± 


14 
04 


130 ±27 ±14 
27 ± 6 ± 3 


5.5 


0.80tg;?o ± °- 03 


8.4 ±1.8 ±1.0 -0.23 ±0.19 ±0.06 




12.2 ±0.7 


1.00 ±0 


06 


39 ± 8 ± 4 






5 = 3.59 ± 0.19 ± 0.12 A<5 = -0.05 ± 0.19 ± 0.06 




-> A' ± 7r + 7r- 


24.7 ± 1.3 


1.68 ±0 


12 


64 ± 14 ± 7 












— > K TV 


10.9 ±0.6 
12.8 ±0.7 


3.33 ±0 
1.24 ±0 
2.09 ±0 


13 
07 
12 


128 ±21 ±12 
48 ± 8 ± 4 
80 ± 13 ±8 


8.2 






8.3 ±1.4 ±0.8 ±0.04 ±0.15 ±0.04 


^(1430)* 
















7.0 ± 1.3 ±0.9 




<^A'*(1410) ± 


28.0 ±2.2 


5.71 ±0 


44 


64 ± 3lia? 


< 2 






< 4.3 (2.4 ± 1.2 




tpK^VHOif 


20.8 ± 1.4 


2.27 ±0 


16 


(90 ±32 tH) 1 


< 2 






< 15.0 




<pK 2 (1820)* 


21.6 ± 1.5 


2.35 ±0 


18 


122 ± 40±^ 


< 2 






< 16.3 





resonance, and find consistent results. The difference is 
accounted as a systematics uncertainty, while the yield of 
decay B ± — > 95i4T2(1770) ± is used to obtain its branching 
fraction. We quote only upper limits on the two branch- 
ing fractions as their correlation is not accounted for in 
the central values. For the tpK 2 and tp K* (1410) decays, 
we vary the longitudinal polarization fraction between 0.5 
and 0.93, and constrain it to 0.8 in the nominal fit. Polar- 
ization variations are included in the branching fraction 
calculations. We vary the kinematic parameter describ- 
ing — ► K ± 7t + tt~ decay (ro2 in Ref. for various 
partial waves of the quasi-two-body K% decay channels 
and take the largest variations as the systematic uncer- 
tainties. The systematic uncertainties in efficiencies are 
dominated by those in particle identification, track find- 
ing, and Kg and ir° selection. Other systematic effects 
arise from event-selection criteria, tp and Kj branching 
fractions, and number of B mesons. 

In summary, we have performed an amplitude analysis 
and searched for CP-violation with the £? — > tpKj 
decays which include significant Ki(1270) and iv|(1430) 
contributions. Our results are summarized in Table [U 
The polarization measurement in the vector-tensor B^ 
decay is consistent with our earlier measurement in the 
B a — > tp if2(1430) decay [l[ and with the naive expec- 



tation of the longitudinal polarization dominance. How- 
ever, our first measurement of polarization in a vector- 
axial-vector B meson decay indicates a large fraction of 
transverse amplitude, similar to polarization observed in 
the vector-vector final state B -> tpK* (892) [HQ. Both 
measurements indicate substantial At+i (or still possible 
A1-1 for vector-axial-vector decay) amplitude from an 
uncertain source. Among potential sources are penguin 
annihilation, electroweak penguin, QCD rescattering, or 
physics beyond the standard model [3j]. 
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FIG. 2: Left column: projections onto the variables rriK-mv 
(a), AE (b), Hi (c), and H 2 (d) for the signal <^A'i(1270) ± 
candidate. Right column: projections onto the variables rriKw 
(e), AE (f), Hi (g), and H 2 (h) for the signal ^ (1430) ± 
and ip(K7r)g candidates combined. The step in (g) is due 
to selection requirement Hi < 0.6 in the channel with n°. 
Data distributions are shown with a requirement on the 
signal-to-background probability ratio calculated with the 
plotted variable excluded. The solid (dotted) lines show the 
signal-plus-background (combinatorial background) PDF pro- 
jections, while the dashed lines show the full PDF projections 
excluding ipKf (left) or ipK^lAm)* 1 (right). 
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